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Abstract First data on inclusive particle production measured in proton-proton collisions at the Large Hadron 
Collider (LHC) are compared to predictions of various hadron-interaction Monte Carlos (QGSJET, EPOS and 
SIBYLL) used commonly in high-energy cosmic-ray physics. While reasonable overall agreement is found for 
some of the models, none of them reproduces consistently the y/s evolution of all the measured observables. 
We discuss the implications of the new LHC data for the modeling of the non-perturbative and semihard parton 
dynamics in hadron-hadron and cosmic-rays interactions at the highest energies studied today. 
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1 Introduction 

The highest energy hadronic interactions measured on Earth result from the collision of cosmic rays (CR) - 
protons and nuclei accelerated in various astrophysical sources that propagate through the universe up to the 
so-called Greisen-Zatsepin-Kuzmin (GZK) cutoff around 10^*^ eV |T]|2l - with air nuclei in the upper atmo- 
sphere |3|. The determination of the primary energy and identity (mass) of such ultra-high-energy cosmic rays 
relies on the study of the cascade of secondary particles, called extensive air-showers (EAS), that they produce in 
the atmosphere |4 |, and its comparison to simulations that include the modeling of hadronic interactions at cm. 
energies {^/s^ w 400 TeV) more than two orders of magnitude higher than those studied at particle colliders 
before the LHC. The dominant source of uncertainty in the interpretation of the highest-energy EAS data stems 
from our limitations to model particle production in strongly interacting systems. Indeed, even at asymptotically 
large energies the collision between two hadronic objects is sensitive to non-perturbative - hadronization, beam 
remnants, soft "peripheral" diffractive scatterings - or semi-hard - saturation of gluon densities, multi-parton 
interactions - dynamics that need still to be constrained directly from experimental data. 

The first LHC data have extended by more than a factor of three the cm. energies for which we have direct 
proton-proton measurements available to test and constrain the ingredients of the hadronic Monte Carlo (MC) 
codes used in CR physics. In this work we compare the predictions of several CR MCs with various inclusive 
observables measured at the LHC which are sensitive to non-perturbative and semihard QCD dynamics: 

(i) total inelastic p-p cross sections CT/,,^/, 

(ii) pseudorapidity density of charged particles at midrapidity dNch/dr]\r]=Q, 

(iii) event-by-event distribution of the charged particle multiplicity P(Nch), 
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(iv) energy distribution of (very) forward particles JA'/JZ?y||,j|>io.94 and dEi,ad/dri\^^^^T,_^, 

(v) average transverse momentum of the produced hadrons {p±). 

The implications of the data-theory comparisons for the improvement of the description of multiparticle produc- 
tion in the hadronic event generators and for the interpretation of CR results are discussed. The interested reader 
can find more details in 0|6l. 

2 Hadronic collisions at multi-TeV energies 

The inclusive production of particles in high-energy hadronic collisions receives contributions from "soft" and 
"hard" interactions between the partonic constituents of the colliding hadrons. Soft (resp. hard) processes involve 
mainly f-channel partons of virtualities typically below (resp. above) a scale of a few GeV^. 

Soft scatterings give rise to production of hadrons with small transverse momenta and dominate hadronic col- 
lisions at low energies (y^ < 20 GeV). Although soft processes have a virtuality scale not far from Aqcd ~ 0.2 GeV 
and thus cannot be treated within perturbative QCD (pQCD), predictions based on basic quantum field-theory 
principles - such as unitarity and analyticity of scattering amplitudes - as implemented in the Gribov's Reggeon 
Field Theory (RFT) \T\, give a decent account of their cross sections in terms of the exchange of virtual quasi- 
particle states (Pomerons and Reggeons). At high energies the dominant soft contributions are from diffractive 
scatterings where one or both colliding hadrons survive the interaction and few particles are produced. 

(Semi)hard parton-parton scatterings dominate the inelastic hadron production cross-sections for cm. energies 
above a few hundreds of GeV. Hard processes with large |?^| 2> Ag^^, can be treated within perturbative QCD 
in a coUinear-factorized approach in terms of parton distribution functions (PDFs) in the hadron convoluted with 
the elementary parton-parton subprocess computable at a given order in the strong coupling constant as{q^). 
The scattered quarks and gluons produce then coUimated bunches of final-state hadrons (jets) in a branching 
process dominated by perturbative parton splittings described by the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi 
(DGLAP) equations lf8]|9l llOI . and by non-perturbative hadronization (e.g. based on the Lund string model 1111 ) 
when the parton virtuality is below GeV). At increasingly larger cm. energies, one needs to account for 
multi-parton scatterings and parton saturation effects. First, the cross section predicted by the (semi)hard pro- 
cesses exceeds the total inelastic p-p cross section for values of a few GeV indicating that multiple parton 
interactions (MPI) occur per collision. Second, for decreasing but still perturbative values, parton scatterings 
receive major contributions from the region of low fractional momenta (x = /7,„,„„„/p,,„rf,„„), where the gluon distri- 
bution rises very fast. In this regime, around a "saturation scale" Q^,,, of a few GeV^, parton branching and fusion 
processes should start to compensate each other saturating the growth of the PDFs as x — >■ Ii l2l . 

The MC event generators of high-energy hadronic collisions used in CR physics - such as QGSJEXOI and 
n I13II14I . SIBYLL lll5lll6lfT7ll and EPOS |18| - have evolved starting up from the RFT approach, based on 
Pomeron degrees of freedom and thus naturally accounting for soft dynamics, generalized to include perturba- 
tive parton-parton processes via "cut (hard) Pomerons" diagrams. Multi-scattering phenomena (gluon saturation, 
MPI) are also implemented through various procedures [SJ. 

3 LHC data versus cosmic ray MCs 

Extensive air showers initiated by interactions of primary cosmic-ray particles with air nuclei in the atmosphere 
constitute multi-step cascade processes involving electromagnetic and hadronic processes. The electromagnetic 
part is well described theoretically, whereas the hadronic interactions are modeled as summarized in the previous 
Section. The most important EAS observables are the the depth in the atmosphere where the number of charged 
particles reaches its maximum Xniax(g/cm^), the number of particles at maximum A'max, and the number of elec- 
tromagnetic particles (e*, y) and muons (^X*) at ground |T9l . 

The relation between high-energy hadronic interactions and EAS observables has been studied numerically 
in detail in [20]. The depth of shower maximum X^^x depends mainly on (i) the inelastic cross section {Oinei) of 
the primary particle with air nuclei, (ii) the corresponding energy fraction (inelasticity) transferred to secondary 
particles but the most energetic "leading" one emitted at very forward rapidities, relative to the primary particle, 
and (iii) the multiplicity (Nd,) of the primary and subsequent very high-energy interactions, which defines how the 
energy is distributed to secondary particles and corresponding sub-showers (and results in a given A'max)- The rest 
of EAS properties at ground are closely related to X^ax and A'max- We discuss below how the LHC measurements 
constrain the collision-energy evolution of quantities such as (7,>,e;, N^h, or the inelasticity. 
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3.1 Inelastic p-p cross section 



A fundamental quantity of all CR models is the total hadronic cross section Oioi and its separation into elastic 
and inelastic (and, in particular, diffractive) components. The measurement of Ogi is accessible thanks to various 
forward proton detectors such as TOTEM 1 21 1 and ALFA (ATLAS) fT2\ in the LHC tunnel area. The ATLAS f23l 
and CMS 1 24 1 experiments have already reported a value Oinei ~ 70 mb at ^/s = 7 TeV which seems to favour the 
lowest of the two (inconsistent) values previously measured at Tevatron (1.8 TeV) pointing to a slightly slower 
increase of the hadron production cross sections with cm. energy as in EPOS and QGSJEtOI (Fig.[T]left). 




Fig. 1 Collision-energy dependence of the inelastic p-pip-p) cross section (left) and of the midrapidity charged hadron multiplicity 
density (right) including the latest LHC data and cosmic-ray MCs predictions. 



3.2 Charged particle multiplicity 



The charged hadron pseudorapidity density at LHC energies provides an important constraint for the modeling 
of the redistribution of energy in the first CR interactions in the atmosphere. The midrapidity measurements of 
dNch/drj |,,=o « 3.5, 4.5, 6at^ = 0.9, 2.36, 7 TeV by ALICE [251, ATLAS \26\ and CMS |271 indicate that the 
multiplicity changes smoothly in the lab energy range from 4 x lO'^ to 3 x lO'* eV, being well reproduced (with 
the exception of EPOS) by the current interaction models used for EAS simulations (Fig. [T] right). In addition, 
first results on the heavy-ion (Pb-Pb) multiplicity and its centrality dependence at ^/s^ = 2.76 TeV II28I provide 
extra important cross-checks on the role of initial-state gluon saturation effects in collisions involving nuclei, 
such as those of CR with air. 



3.3 Multiplicity probability distributions 



The multiplicity distribution P{Nch), i.e. the probability to produce Nch charged hadrons in an event provides 
important differential constraints on the internal details of the hadronic interaction models. The low multiplicity 
part is mostly dominated by the contributions from diffraction (single-cut Pomeron exchanges), whereas the 
tail of the distribution gives information on the relative contribution of multiparton scatterings (multi-Pomeron 
exchanges). The experimental measurements (Fig. |2]i at the three cm. energies measured so far at the LHC, 
indicate that the high-Nd, tail (left) is underestimated by EPOS and QGSJETOI, whereas SIBYLL and QGSJETlI 
get a bit closer, sometimes overestimating the data. In the low-Nch region around P(Nch) ~ 4 (right panel), only 
EPOS globally reproduces the experimental results whereas the rest of the models overestimate the measurements 
up to +30% for SIBYLL. The peak is even shifted towards lower multiplicity in the case of both QGSJET models. 
Thus, even if the average p-p multiplicities are well reproduced by most models (Fig. [T] right), the details of 
their probability distributions are missed and indicate possible paths for improvement of the different model 
ingredients. 
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Fig. 2 Multiplicity distributions of charged hadrons, P(Nch ) , measured by ALICE in p-p events at y'J = 0.9, 2.36 and 7 TeV ( 25 1 (left) 
compared to QGSJETOl and II, SIBYLL, and EPOS. The right plot shows a "zoom" in the low multiplicity range for i/J = 7 TeV. 



3.4 Forward energy flow 

A key quantity for air shower development is the distribution of neutral energy emitted at very forward rapidities 
as it provides constraints on the production of leading hadrons (inelasticity) as well as on the transfer of energy 
from the hadronic core to the electromagnetic cascade (via — >■ yy decays). The recent LHCf measurement 
of the 7 spectrum for rapidities above |r)| w 8.8 |29| is compared to model simulations in Fig.[3](left) |6|. The 
simulations are in relatively good agreement with the data within the systematical uncertainty (not shown here) 
although for Ey< 1.5 TeV the spectrum slope is harder in the data than in the predictions. Of similar importance is 
the measurement of the energy flow and particle spectra in the forward range | r] | = 5 - 10 fSO]. This is an angular 
range that has been historically very difficult to access in collider experiments but that is partially covered by 
various detectors at the LHC such as TOTEM (2D and CASTOR (CMS) (H]. Preliminary results of the energy 
flow in the \ r] \ =3-5 range indicate a good data-model agreement at -y/y = 0.9 and 7 TeV I32II (Fig.[3]right). 




Energy (GeV) 3 3-2 3.4 3.6 3.8 4 4.2 4.4 4.6 4.8 



Fig. 3 Forward energy distributions in p-p collisions at the LHC compai'ed to CR Monte Carlos, for photons at | ?] | > 10.94 measured 
by LHCf at 7 TeV [251 (left) and for inclusive hadrons at 3 < | J] | < 5 measured by CMS at 0.9, 7 TeV r32 1 (right). 
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3.5 Transverse momenta of hadrons 

Although the transverse momentum spectrum of the produced hadrons (or its average p±) at the LHC does not 
have a direct impact on the interpretation of air shower data - the lateral distributions of particles at ground is 
rather defined by multiple Coulomb scattering and by the spectra of secondaries at much lower energies - such 
a measurement is of importance for checking the overall physics consistency of soft and hard interaction mech- 
anisms implemented in the models. Indeed, at high energies the peak of the perturbative production comes from 
interactions between partons whose transverse momentum is around the saturation scale, p± ~ Q^at , producing 
(mini)jets of a few GeV which fragment into hadrons. In models with saturation of parton densities, the mean 
transverse momentiun of the produced hadrons is of the order of the saturation scale Q^ai in the high-energy limit. 

In Fig. |4] (left) we show the energy evolution of the mean p^ measured experimentally compared to the CR 
event generators and to the PYTHIA 8 MC 133]. All the RFT MCs but EPOS predict a very moderate increase of 
(p±) with energy, reaching (p^) ~ 0.6 GeV/c at GZK energies which is only 0.05 GeV/c above the current CMS 
result at 7 TeV, reflecting the moderate assumptions made on the saturation of the low-x parton densities. EPOS 
predicts a significantly larger {pl) ^.^k ^ ^ GeV/c due to the inclusion of final-state collective parton expansion 
effects. PYTHIA 8 - whose dynamics is dominated by (mini)jet production with a running p±^ cutoff that mimics 
parton saturation effects (5) - predicts a higher average (p^) w 1.2 GeV/c. 
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Fig. 4 Left: Average of chai'ged pailicles at midrapidity in p-p (p-p) collisions as a function of y'i compared to the predictions 
of CR event generators and PYTHIA 8. Right: Collision-energy dependence of the ratio of antiprotons to n'^ yields at midrapidity 
predicted by various CR models. 



The level of (dis)agreement between the inclusive hadron observables measured at the LHC discussed in (5j 
|6j and each one of the four CR hadronic MCs considered, is summarized in Table [T] The event generators give 
an overall decent description of all measurements but model improvements, particularly, those related to the 
treatment of inelastic diffraction and of the parton saturation mechanism, are desirable. 
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Table 1 Level of overall agreement between QGSIETOI, QGSIETII, SIBYLL 2.1 and EPOS 1.99 with inclusive charged hadron results 
measured in collisions at 0.9, 2.36 and 7 TeV: inelastic cross section (T„„;, pseudorapidity densities dN^.iJdr\\j^^Q, multiplicity prob- 
abilities P(Nch) for low and high values of N^i,, and mean transverse momentum (pa_). A tick (/) indicates a reasonable data-model 
agreement within experimental uncertainties, and 1^ (J|) that the MC tends to over (under) estimate the measurements. 
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4 Conclusions 

Event generators used in higii-energy cosmic ray (CR) physics include a description of hadronic interactions that 
partially depends on non- and semi-perturbative QCD dynamics that need to be calibrated with experimental data. 
The highest energy cosmic rays measured on Earth at Eiat, w 10^" GeV collide with the nuclei in the atmosphere 
at cm. energies more than two orders of magnitude above those studied at particle colliders before the LHC. 
The most recent results from p-p collisions at the LHC (equivalent to lab energies around Ei„i, = 3 x 10'* eV) 
are of big help to constrain the details of multiparticle production in the Monte Carlos used to describe CR air 
showers. The measured characteristics of the bulk of hadron production at multi-TeV energies does not reveal 
serious deficiencies in any of the models. This gives a strong support to the interpretation of the results in the CR 
"knee" energy range (£/;,/, = lO'^^eV) in terms of conventional primary spectrum and nuclear mass composition 
and disfavours some proposed speculative ideas that the change of the CR spectral slope could be due to a sudden 
change in the hadronic interaction mechanism above 2TeV cm. energy (see e.g. 1341 "). 

Although the first LHC measurements support a conventional extrapolation of the known features of multi- 
particle production to the highest known energies, none of the models is in perfect agreement with all the hadronic 
observables measured at the LHC (see Table [T}. In particular, extrapolations at the GZK-cutoff energies span a 
range of predictions - e.g. dNch/dr\\^=(, w 10 (EPOS, SIBYLL) - 50 (QGSJETlI) for the particle densities and 
(px) ~ 0.6 (SIBYLL, QGSJEtOI) - 1 (EPOS) GeV/c for the mean hadron transverse momentum - that justify the 
concurrent use of various MCs to gauge the uncertainties connected to hadronic interaction models in the inter- 
pretation of the cosmic ray data. New retuning of model parameters and reconsideration of model assumptions, 
e.g. for EPOS lt35J and QGSJET t36j, are currently underway. 

Further improvement of our understanding of the strong interaction and of the properties of cosmic-rays at 
multi-TeV energies will be provided by the LHC with p-p data at the nominal cm. energy of = 14TeV 
(corresponding to CR protons of lO'^ eV in the lab frame), as well as from the expected proton-nucleus {p- 
Pb) runs at ^/s^ =8.8 TeV |j37j, since the CR-induced air showers in the upper atmosphere are mostly from 
(p,a,Fe)H-(N,0) collisions. Also, since baryon-induced subshowers lead to a higher number of muons at ground 
than meson-induced ones, the energy dependence of the baryon production rate (see e.g. Fig. |4] right), its relation 
to the centrality of the collision, and the momentum distribution of the baryons constitute also important quantities 
to be measured. 
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